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Abstract 

Detailed measurements of the size of roughness ele- 
ments on ice accreted on models in the NASA Lewis 
Icing Research Tunnel (IRT) were made in a previous 
study. Only limited data from that study have been 
published, but included were the roughness element 
height, diameter and spacing. In the present study, the 
height and spacing data were found to correlate with 
the element diameter, and the diameter was found to be 
a function primarily of the non-dimensional parame- 
ters freezing fraction and accumulation parameter. The 
width of the smooth zone which forms at the leading 
edge of the model was found to decrease with increas- 
ing accumulation parameter. Although pr eliminar y, 
the success of these correlations suggests that it may be 
possible to develop simple relationships between ice 
roughness and icing conditions for use in ice-accretion- 
prediction codes. These codes now require an ice- 
roughness estimate to determine convective heat 
transfer. Studies using a 7.6-cm-diameter cylinder and 
a 53.3-cm-chord NACA 0012 airfoil were also per- 
formed in which a Vi-min icing spray at an in itia l set of 
conditions was followed by a 9V£-min spray at a second 
set of conditions. The resulting ice shape was com- 
pared with that from a full 10-min spray at the second 
set of conditions. The initial ice accumulation ap- 
peared to have no effect on the final ice shape. From 
this result, it would appear that accreting ice is affected 
very little by the initial roughness or shape features. 

Nomenclature 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Airfoil leading-edge or cylinder diameter, cm 

c Pta Specific heat of air, cal/g K 

c PtW Specific heat of water, cal/g K 

d Roughness-element diameter, cm 

evap Evaporation terms in energy equation, °C 
h Roughness-element height, cm 

h c Convective heat-transfer coefficient, 

cal/s cm 2 K 


LWC Cloud liquid-water content, g/m 3 

MVD Water droplet median volume diameter, \xm 
n Freezing fraction, dimensionless 

s Spacing between roughness elements, cm 

V Airspeed, m/s 

tf Freezing temperature of water, °C 

Temperature at surface of ice, °C 
tjt Static temperature, °C 

it* Total temperature, °C 

w Smooth-zone width, cm 

fi 0 Stagnation-zone droplet collection efficiency, 

dimensi onless 

8 Droplet median volume diameter, jjm 

^ Droplet energy transfer terms in energy bal- 
ance, °C 

A f Latent heat of freezing of water, cal/g 

pt Ice density, g/m 3 

0 Air energy transfer terms in energy balance, 
°C 


Introduction 

This paper presents the results of a study of the effects 
and characteristics of the roughness which forms on 
the surface of a glaze-ice accretion. The freezing rate 
of water im ping in g to form glaze ice is strongly de- 
pendent on the rate of heat transfer. Because surface 
roughness is known to have an important effect on the 
convective heat transfer, an estimate of the roughness 
size is presently used to establish heat-transfer coeffi- 
cients in ice-accretion prediction codes such as the 
NASA Lewis LEWICE code 1 ' 2 . In this paper, correla- 
tions of measurements of roughness made by Shin 3 will 
be presented. In addition, the effect of initial ice 
roughness and ice formation on the final ice shape was 
explored, and these results will be presented. 

In the past few years a number of studies have been 
mad* to try to understand the physical processes which 
occur during an icing encounter of an aircraft in flight 



This interest in icing physics has been a response to the 
need to improve the predictive capability of analytical 
ice-accretion models such as the NASA Lewis LE- 
WICE 1 - 2 code. 

The first of these studies was made by Olsen and 
Walker 4,5 who studied ice accretion with high-speed 
motion pictures. Based on their observations, they 
questioned the traditional model of ice accretion pro- 
posed by Messinger 6 in 1953. The Messinger model 
for gh« ice assnm^H that water which did not freeze 
on impart ran back along the surface to freeze farther 
hark Olsen and Walker observed no runback of a liq- 
uid film but rather noticed that unfrozen water coa- 
lesced into beads on the surface. These beads subse- 
quently froze into closely spaced, approximately hemi- 
spherical roughness elements. In the early stages of 
icing exposure, a relatively smooth zone at the leading 
edge of the test model was also observed; this zone 
OTtmHwt back some distance until a sudden transition 
to the rough zone occurred. 

Additional close-up studies of ice accretion were later 
performed by Hansman, Reehorst and Sims 7 using 
high-speed video techniques. They noted the forma- 
tion of a thin rime ice layer on the surface of the model 
when icing started. This phenomenon was attributed 
to conduction heat transfer into the model, so its occur- 
rence could be dependent on model construction. The 
authors pointed out that once such a layer had formed, 
any subsequent ice formation would be unaffected by 
the model surface treatment (polishing, waxing, etc.) 
that might be done to influenc e the initial bead forma- 
tion of water. The growth of water beads through coa- 
lescence was noted to occur until the bead was suffi- 
ciently large to be influenced by aerodynamic farces. 
When this occurred, the bead was suddenly swept 
downstream where it froze. 

The subsequent freezing observed in reference 7 was 
attributed to a higher heat transfer rate where the 
boundary layer was tripped. However, Bragg and co- 
workers 8,9 found that roughness elements of the size 
measured by Shin 3 did not cause immediate transition 
of the boundary layer from laminar to turbulent In- 
stead, while the transition began at the edge of the 
rough zone, si gnificant d j stance along the surface was 
required before fully turbulent flow was established 
well bade on the airfoil. 

The formation of the hemispherical beads of roughness 
through coalescence of liquid water on the surface of 
the ice would be strongly influenced by water surface 
tension Hansman and Tumock 10,11 studied bead for- 


mation with water and with a water-surfactant mixture. 
The addition of surfactant to the spray water reduced 
the surface tension of the water by a factor of 2 and 
produced smaller beads. The final ice formation which 
resulted with surfactant addition was opaque and had 
more distinct horns than that from untreated water. 
S imilar results with surfactant were reported by oth- 
ers. 12 These studies suggested that surface effects, in 
addition to convective heat transfer, may have a strong 
influence on the final ice shape Surface phenomena 
are outside the scope of this paper, but deserve more 
investigation. 

The LEWICE 2 ice-accretion code currently es timat e s 
the roughness element size by balancing aerodynamic 
forces, which tend to move the droplet along the sur- 
face, with surface-tension forces, which tend to restrain 
the droplet motion. The roughness size is then used to 
determine the convective heat transfer. Shin’s 3 objec- 
tive was to establish relationships between icing condi- 
tions and the roughness characteristics. If the rough- 
ness correlated with icing conditions in a simple way, 
the roughness size could be calculated from the corre- 
lations in the LEWICE code. Because the present pa- 
per draws heavily on the work of Shin, that research 
and its results will be discussed in the next section in 
some detail. 

It bag been speculated that the initial toughness which 
forms during ice accretion may have a strong effect on 
the final ice shaqie. Roughness features could provide 
locally enhanced collection efficiencies for water 
droplet impingement, and they would also tend to in- 
crease the local heat transfer for higher freezing rates 
In this way, small ro ughness features would grow rap- 
idly, they would dominate the accretion process and 
deter mine the final ice shape. Hangman and Tur- 
nock," however, found that, “as the accretion grows, 
the effect of initial conditions begins to wash out,” and 
the initial surface condition had no effect on final ice 
shape. The second objective of the present study, then, 
was to explore further the effect of initial roughness on 
final ice shape. 

The present study used a 7.6-cm- (3-in-) diameter hol- 
low al uminum cylinder and a 53. 3 -cm- (21-in-) chord 
solid al uminum NACA 0012 airfoil to study the effect 
of initial spray conditions on final ice shape. Ice 
shapes were recorded for 10-min sprays in which the 
initial 14 min was at a different set of conditions from 
the final 9 'A. These shapes were compared with those 
obtained for a full 10 -min spray at the final set of con- 
ditions. Icing conditions included static temperatures 
of -6°C to -1°C (21°F to 30°F), airspeeds of 67 and 89 
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m/s (150 and 200 mph), a median volume diameter 
drop size of 20 pm and liquid-water contents of .5 and 
1 g/m 3 . Tests were performed in the NASA Lewis Ic- 
ing Research Tunnel (IRT). 


Correlations from Shin’s 3 St 


Shin used both a digital and 35-mm camera to photo- 
graph the roughness characteristics on the leading edge 
of a 53.3-cm- (21-in-) chord NACA 0012 airfoil after 
s pray s of 1 to 3 min utes in conditions which would 
produce glaze ice. All tests were mad e in the IRT. 


Figure 1 is a sketch which shows features of a typical 
glaze-ice accretion beginning to form after a short ex- 
posure to an icing spray. From digital photographs of 
such ice accretions for a variety of test conditions and 
exposure times Shin found that the ice-roughness ele- 
ments took the form of hemispheres on the surface of 
an ice substrate. These roughness elements formed 
some distance bade from the leading edge leaving a 
smooth zone whose width decreased with exposure 
♦imp Shin measured the diameter, d, height, h, and 
s paring c tn characterize the element roughness and 
recorded the width, w, of the smooth zone. 


Figure 2 shows how the roughness-element spacing, s, 
and height, h, correlate with the diameter, d. The ele- 



gized the importance of developing a good under- 
standing of how roughness forms. 

Shin Sho wed that ice roughness dement size grew with 
time up to 2 min then became constant. A reasonable 
pnctiiiation is that the size of the coalesced beads is 
Hrtrrminrd by the rates of water impingement and 
freezing. For a given rate of water accumulation, a 
high rate of freezing would halt the growth of beads at 
a smaller size than would occur with a low rate of 
freezing. With low freezing rates, the bead size may be 
limitM hy the houndarv-laver thickness: Beads which 
grow sufficiently to protrude from the boundary layer 
are subj ect ed to aerodynamic forces which would tend 
to sweep them downstream, as observed in the studies 
of references 4,5 and 7. In the present work, no at- 
tempt was made to relate roughness size with bound- 
ary-layer thickness; instead. Shin’s roughness dimen- 
sions will be shown here to be correlated with the non- 
dimensional terms familiar to icing research, accumu- 
lation parameter and freezing fraction. 

The accumulation parame ter is a non-dimensional pa- 
rameter which is proportional to the rate at which ice 
accretes. Itisgivenby 
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The freezing fraction was defined by Messinger 6 as the 
fraction of water which froze in the region of im- 
pingement. Its value can be found from an energy bal- 
ance made on the surface at the leading edge of a 
model exposed to an icing cloud: 



Figure 1. Roughness Features on Leading Edge of 
Airfoil. 

ments almost touched so that the spacing was slightly 
greater than one diameter. The height was approxi- 
mately one-half the diameter. Because of the strong 
correlation of spacing and height with the diameter 
shown in Figure 2, it should be adequate to use simply 
the Hiamrter to describe the effects of icing conditions 
on the roughness characteristics. The roughness char- 
acteristics measured by Shin often differed signifi- 
cantly from those predicted by the then-current version 
of the LEWICE code; and this result further empha- 


where ft is a collection of terms representing the drop- 
let energy transfer 

, v ! 

*** p,w 

and 0 represents the air energy transfer: 

V 2 

0 = t„ rf -t a -r— + evap 

^ P * 
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where evap represents a water evaporation term, which 
is usually small. Both ^ and 0 have units of tempera- 
ture. b is a non- dimensi onal term called the relative 
heat factor, it is the ratio of the heat carried to the sur- 
face by im ping in g water to that convected away by the 
airstream. It is defined as: 

, LWCVP»c r . 

K 

The effect of accumulation parameter on roughness 
size is given in figure 3. For data with a fr ee zin g frac- 
tion of approximately .3, Shin’s data in d ica te d that the 
roughness size increased linearly with accumulation 
parameter until A c reached a value of. 1; iorA c greater 
than . 1 , element size was nearly constant 

In figure 4 roughness-element diameters resulting from 
tests with A c > . 1 (i.e., diameter independent of A c ) 
have been plotted as a function of the freezing fraction. 
Roughness size decreased with increasing freezing 
fraction, with the rate of decrease greatest at the lowest 
freezing fractions. 

Finally, Shin’s smooth-zone width has been correlated 
with the accumulation parameter in Figure 5. For 
these conditions, the smooth-zone width was found to 
decrease linearly with accumulation parameter. The 
mechanism by which the smooth zone contracts has 
not yet been fully explained. The data in Figure 5 were 
for a freezing fraction of .3; insufficient data exists to 
determine if other values of freezing fraction would 
give different correlations. 

Description of Experiments 

NASA Lewis Icing Research Tunnel The IRT is 
shown in Figure 6. It has been described in reference 
13. The IRT has a test section width of 2.7 m (9 ft) 
and a height of 1.8m(6ft) It is capable of operation 
at test-section airspeeds up to 160 m/s (350 mph) with 
a blockage of 5% in the test section. A refrigeration 
system permits accurate control of the test-section tem- 
perature from -30 to 1°C (-20 to 33°F) within ± .5°C 
(± 1°F). A water-spray system with 8 spray bars has 
been calibrated to provide controlled test-section liq- 
uid-water contents from .2 to 3 g/m^ and droplet me- 
dian vol ume diameters from 15 to 40 fim. 

Two sets of spray nozzles, known as the mod-1 and 
standard nozzles, are used in the IRT to provide differ- 
ent ranges of liquid-water content and droplet size 13 . 


For the tests of this study, only the mod-1 nozzles were 
used. 

Test Models Two models were used to study the effect 
of initial ice roughness on final ice shape. The first 
was a 7.6-cm- (3-in-) diameter hollow aluminum cyl- 
inder. The cylinder was mounted on stands in the 
center of the IRT test section as shown in figure 7. To 
permit the IRT spray-bar conditions to reach steady- 
state before exposing the model to the icing cloud and, 
thus, to control the exposure time precisely, a shield 
was lowered in front of the model as shown in the fig- 
ure. Once the spray had stabilized, the shield was 
raised; the exposure time was counted from the raising 
of the shield to the completion of the spray. 

The second model which was used to look at the effect 
of initial rou ghness was a 53.3-cm (21-in) NACA 0012 
airfoil mounted vertically in the center of the IRT test 
section. It was a 1.8-m- (6-ft-) span model machined 
from solid aluminum. 

Figure 8(a) shows the airfoil installed in the test sec- 
tion. The model was mounted to the IRT turntable at 
its base and to the ceiling with a pivot Although this 
mounting gave the capability of varying angle of at- 
tack, all tests were run at 0°. The metal shield de- 
scribed above for the cylinder tests could not be used 
with this airfoil because of mounting constraints. In- 
stead, a fabric shield was attached to the leading edge 
with velcro fasteners as shown in Figure 8(b). A rope 
att ache to the fabric was passed through a hole in the 
tunnel ceiling. When spray conditions had stabilized, 
the rope was manually jerked to pull the shield clear of 
the model. An advantage of this method of protecting 
the airfoil during spray stabilization was that the test- 
section blockage was unchanged when the protection 
was removed; consequently, the test-section velocity 
remained constant 

Total Temperature Measurement The IRT total tem- 
perature is measured as the average of 11 type-T ther- 
mocouples distributed over the downstream plane of 
the corner-D t ur nin g vanes. Corner D is just upstream 
of the spray bars. In addition to these measurements, a 
heated Rosemount total -temperature probe was 
mounted on a horizontally traversing mechanism 
mounted just downstream of the test section for the 
airfoil tests. This probe was located to sense the tem- 
perature mid way between the floor and ceiling of the 
tunnel. The probe was traversed across the tunnel be- 
fore and after selected tests to calibrate the comer-D 
average temperatures with the actual test-section tem- 
perature. The traversing temperature probe was found 
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to be mngi stent with the comer-D average temperature. 
For the airfoil tests only, the comer-D temperatures 

were recorded every 2 sec throughout each test Aver- 
age run temperatures were then computed. These av- 
erages agreed with the temperature set by the tunnel 
operators within ±_3°C (.5^). Because of this close 
agreement, the test temperature reported here is the 
temperature requested of the tunnel operators. 

The thermocouple readings were considered accurate 
to within ±.5°C (1°F). With possible additional dis- 
crepancy between the requested and actual average run 
temperatures of ±.3°C ( 5°F), the total uncertainty in 
temperature was about ±.8°C (1.4°F). 

Test-Section Velocity Measurement The IRT uses two 
pitot-static probes, one on each side of the entrance to 
tlv» test section, to determine velocity. For the airfoil 
tests, the total and static pressures from each probe 
were recorded at 2-sec intervals during each test The 
test-section velocity was calculated from the pressures 
for each probe at each recording interval. The run- 
average velocity for each probe was then obtained, and 
these two velocities were averaged to give the test- 
section fiee-stream velocity. The individual averages 
found from the two probes differed consistently by 
about 1 m/s (2 mph). Considering the accuracy of the 
transducer s along with this uncertainty gives an esti- 
mated uncertainty in the velocity measurement of less 
than ±4%. The measured average test section velocity 
never differed from the requested set point by more 
than .5 m/s (1 mph); thus, the velocities reported will 
be the set conditions. 

Sorav-Bar Conditions The IRT is routinely calibrated 
on a regular basis. The cloud median volume diameter 
(MVD) and liquid-water content ( LWC ) were deter- 
mined by applying this calibration to measured spray- 
bar water and air pressures. For the airfoil tests, pres- 
sures were recorded at 2-sec intervals, the MVD and 
LWC calculated from the pressures at each interval 
and, finally , the average run MVD and LWC computed. 
These average run values agreed with the requested 
values typically within 1%, so the requested values will 
be reported here for all tests. Considering both the in- 
accuracy of the transducers used in the measurements 
and the fluctuations with time, the typical uncertainty 
in the liquid-water content was found to be less than ± 
2.2% and in the median volume diameter, less than ± 
12.7%. 

Test Procedure Tests were performed by first estab- 
lishing the desired test-section velocity and tempera- 
ture. For the cylinder tests, the remotely-operated 


shield was lowered in front of the model at this point; 
the airfoil shield was manually placed on the model 
before starting the tunnel fan. Water spray was then 
initiated The shield was raised when the spray condi- 
tions had stahiii-zed and spray timing was started. 

To deter mine the effect of initial roughness on the final 
ire shape, a series erf paired tests were performed. For 
the first test of the pair, an initial 'A-min accretion was 
allowed to form on the model followed by a 9'/2-min 
accretion with conditions sufficiently different from the 
initial to give different freezing fractions and, thus, 
roughness sizes. The second test of the pair used the 
final conditions from the first test for the full 10-min 
spray. In effect, then, the initial %-min roughness for 
the two tests should have been different The fina l 10- 
min ice shapes for each test were compared. 

The tests were performed by shielding the model after 
the initial Yr-min spray of the first test of the pair. The 
shield was dropped in front of the model remotely for 
the cylinder, but for the airfoil tests, the tunnel was 
brought to idle to permit personnel entry to the test 
section to place the shield on the model. The second 
set of tunnel conditions was then established, and, once 
the new conditions bad stabilized, the shield was lifted 
to mqmge the model to the second set of conditions for 
the aHHitinnal 9‘A min. The tunnel was brought to idle, 
the final ice shape was recorded and the model 
cleaned. For the second test run of the pair, the second 
set of tunnel condi tions was used without interruption 
for a full 10-min icing exposure and the ice shape 
again recorded. After the ice shape was recorded, the 
model was cleaned and the procedure repeated for the 
next spray series. 

The ice shape was recorded by melting a thin gap 
through the ice normal to the model span-wise axis. 
The shape was then traced onto a cardboard template. 
All ice shapes presented here were taken midway be- 
tween the floor and ceiling of the tunnel. These shapes 
were i ^ter digitized for computer storage of the coordi- 
nates for preparation of comparison plots. 

F.ffect of Initial Accretion on Final Ice Shape 

In this section, the ice shapes will be compared for the 

with alternate initial spray conditions. 

fSdinder Tests The first series of tests were performed 
ncing the 7.6-cm-diameter cylinder and metal shield 
described above. 
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In figure 9 the effect of changing the liquid-water con- 
tent is shown. All accretions for the tests of Figure 9 
were nyiria with a static temperature of -4°C (24°F), a 
velocity of 67 m/s (ISO mph) and a median volume 
drop diameter of 20 jim. In 9(a) an ice shape produced 
by graying with a liquid-water content of 1 g/m 3 for l A 
min followed tty .5 g/m 3 for 9!A min is compared with 
two shapes with LWC of .5 g/m 3 for 10 min. Compari- 
son of the latter two ice shapes gives some indication 
of the repeata bility of shapes when the same conditions 
are specif ie d in the IRT. Within normal repeatability 
of ice shapes, there appeared to be no effect of the ini- 
tial V£-min spray on the final shape. In 9(b), the initial 
%-min spray had an LWC of 1 g/m 3 , followed by .5 
g/m 3 for 9% min. This shape is compared in the figure 
with (me produced with a 10-min spray at 1 g/m 3 . 
Again, there was no significant difference in ice shape. 

The freezing fraction was estimated to be .20 for the 
tests with a liquid-water content of 1 g/m 3 , and .35 at 
.5 g/m 3 . From Figure 4, these freezing fractions sug- 
gest that the LWC of 1 g/m 3 would produce a roughness 
about 25% higher than the .5-g/m 3 condition. A com- 
parison of the ice shapes produced by the .5- and 1- 
g/m 3 conditions in Figures 9 (a) and (b) shows a sig- 
nificant difference; yet, in each case, the initial Varvoin 
spray condition had no effect on the final ice shape. 

The effect of icing with two temperatures is illustrated 
in figure 10. The accretions were made with a velocity 
of 67 m/s (150 mph), a median volume drop diameter 
of 20 pm and a liquid-water content of .5 g/m 3 . The 
first ice shape in figure 10(a) had an init i a l l /4-min 
spray at -6°C (21°F) static temperature {n = .51) fol- 
lowed Ity 9 l A min at -4°C (24°F) ( n = .35). This shape 
is compared with two made with 10-min sprays at a 
constant -4°C. The latter two shapes indicated that the 
ability to repeat ice shapes at this condition was quite 
good. The effect of the -6°C initial spray appeared to 
be iti signific ant in determining the final ice shape. 

Figure 10(b) shows ice shapes accreted for the same 
conditions of 10(a) except that the spray-condition se- 
quence was reversed, with the first shape made with a 
spray at -4°C (24°F) for l A min followed by -6°C 
(21°F) for 9 X A min. This shape only approximates that 
from a 10-min spray at -6°C. The reason for this poor 
nrntrh is not apparent but will be discussed further af- 
ter presenting additional cylinder test results. 

In the next series of tests, static temperatures of -6 and 
-3°C (21 and 26°F) were used. Test results are re- 
ported in Figures 10(c) and 10(d) for sprays with the 
sjime velocity, MVD 9 and LWC as used in Figures 10(a) 


and (b). In 10(c), the first shape resulted from an ini- 
tial !4-min spray at a static temperature of -6°C (n = 

.51) followed Ity 9 l A min at -3°C (n = .24). The com- 
parison with a shape produced with a 10-min spray at - 
3°C is again very good Reasonable, although not 
perfect, agreement was found between an ice shape 
generated with an initial spray at -3°C for l A min then - 
6°C for an additional 9 l A min and a shape from a 10- 
min spray at -6°C (figure 10(d).) 

Both Figures 10(b) and 10(d) represent cases in which 
a low-freezing-fraction spray was followed by a higher 
one. In each case, the ice-shape did not fully match 
that for a 10-min spray at the higher freezing fraction. 
This result might suggest that the larger roughness 
elements produced at the lower freezing fraction had a 
discernible effect on the final shape. However, Figure 
9(a) showed that a l A-min spray with a freezing frac- 
tion of .2 followed by 9 l A min at n = .35 produced a 
shape which was an excellent match for a 10-min 
shape at n = .35. Furthermore, the shapes of Figure 
10(d) are a better match than those of 10(b) even 
though the freezing-fraction difference for the 10(d) 
tests was greater. Thus, further studies are needed to 
determine if foe umnatriiftri ice shapes seen in Figures 
10(b) and (d) are repeatable, and, if so, to establish a 
physical explanation. 

Airfoil Tests The second group of tests used the 53.3- 
cm- (21-in-) chord NACA 0012 airfoil with the fabric 
shield. The test procedure was the same as for the 
cylinder tests, except th at the fabric shield had to be 
manually placed over foe %-min ice formation at foe 
end of the initial spray. Shield removal involved jerk- 
ing it away using an attached rope as described above 
in the Description of Experiments section. This action 
could potentially have damaged the initial ice forma- 
tion as the fabric tended to drag over the surface. To 
evaluate foe possible extent of disturbance, a few trials 
were madfc in which the ice was observed before and 
after shield removal. No visual effect was apparent 

The ice-shape comparisons for the airfoil tests are 
given in Figure 11. All tests reported in this figure 
usedanA/FDof20 pmand anZ-lFCof .5 g/m 3 . Figure 
11(a) compares a 10-min ice shape made at a tem- 
perature of -6°C (21°F) with the shape resulting from 
first spraying at -3°C (26°F) for Vi min then spraying 
at -6°C (21°F) for 9'A min All temperatures quoted 
are the static temperatures, and the velocity was 67 m/s 
(150 mph) for all the sprays of Figure 11(a). The 
freezing fraction for these conditions was es tim a t ed to 
be .58 at -6°C (21°F) and .27 at -3°C (26°F). Reference 
to Figure 4 suggests that these freezing fractions could 
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have given a 30 - 40% difference in surface roughness. 
Here, again, in contrast with the results of Figures 
10(b) and (d) the test with the high initial roughness 
gives a shape which matches a test with a low initial 
roughness. 

Figure 11(b) looks at the effect of changing the test- 
section velocity during the spray. For these tests, the 
static temperature was maintained at -6°C (21°F). The 
initial i4-min spray was made with a test-section ve- 
locity of 67 m/s (150 mph). This was followed by a 
9*>4-min spray at 89 m/s (200 mph). The final ice 
shape from these sprays is compared with a 10-min 
spray at 89 m/s (200 mph). The lower velocity gave a 
freezing fraction of about .58 at these conditions, and 
the higher velocity, n = .40. The correlation of Figure 
4 suggests that roughness size may differ by about 
30%. Again, no effect of initial spray conditions on 
final ice shape was evident 

Concluding Remarks 

In general, this study showed that 10-min ice shapes 
appear to be independent of initial 14-min icing condi- 
tions. This finding is consistent with the observations 
of Hansman and Turnock, 11 who found that ice shapes 
accreted on surfaces of different materials were the 
same. The initial ice features are covered by subse- 
quent accretion and do not appear to influence later 
accretions significantly. The final ice shape appeared 
to be determined by the spray conditions at which the 
dominant accumulation occurred. 

Using limited data from Shin 3 , this paper successfully 
correlated roughness size with freezing fraction and 
accumulation parameter. If such correlations are pos- 
sible, they can be readily incorporated into ice- 
accretion-prediction codes to permit accurate convec- 
tive heat-transfer coefficient determination for any 
specified icing condition. At present, such codes must 
rely on estimates of roughness size to determine heat 
transfer. The correlations presented here are based on 
limited data and are preliminary. Additional rough- 
ness data are needed for confirmation. 
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Test Section: Room 


1.8m high, 2.7m wide, 1 I — LiUl 

160m/s Velocity 

Figure 6. NASA Lewis Icing Research Tunnel (IRT). 



Figure 7. Test Cylinder and Shield Mounted in IRT. 
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(a) Airfoil Installed in IRT Test Section. (b) Fabric Shield in Place on Leading Edge of Airfoil. 

Figure 8. 53.3-cm-chord NACA 0012 Airfoil. 



Spray 

(a) Initial (1/2 min): LWC = 1 g/m\ n = .20 

Final (9 1/2 vam)-.LWC = .5 g/m 3 , n = .35 



Spray 

Initial (1/2 min): LWC = .5 g/m\ n = .35 
Final (9 1/2 min): LWC = 1 g/m J , n = .20 

LWC = 1 g/m 3 , n = .20 


Full (10 min): LWC = .5 g/m 3 , n = .35 Full (10 min): 

Full (10 min): LWC = .5 g/m 3 , n = .35(Repeal.) 

Figure 9. Effect of Initial LWC on Final Ice Shape for Cylinder. Cylinder Diameter, 7.6 an (3 in). 
V= 61 m/s (1 50 mph), t„ = -4°C (24°F), MVD = 20 pm. 


(b) 
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(a) Initial (1/2 min): /* = -6° C (21°F), n = .51 0>) Initial (1/2 min): /* = -4°C (24°F), n = .35 

Final (9 1/2 min):/* = -4°C (24°F), n = .35 Final (9 1/2 min): /* = -6°C (21°F), n = .51 

Full (10 min): /* = -4°C (24°F), n = .35 Full (10 min): /* = -6°C (21 °F), « = .5 1 

Full (10 min): - -4°C (24°F), n = .35 (Repeat.) 






Spray Spray 

(°) Initial (1/2 min): t tt = -6°C (21 °F), n = .51 W Initial (1/2 min): /* = -3°C (26°F), n = .24 

Final (9 1/2 min): /* = -3°C (26°F), n = .24 Final (9 1/2 min): /* = -6°C (21 °F), n = .51 

Full (10 min): /*= -3°C (26°F), w = .24 Full (10 min): /*= -6°C (21°F), n = .51 

Figure 10. Effect of Initial Temperature on Final Ice Shape for Cylinder. Cylinder Diameter, 7.6 cm (3 in). 

V= 67 m/s (1 50 mph), MVD = 20 pm, LWC = .5 g/in 3 . 
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Initial (1/2 min): t M = -3°C (26°F); n = .27 

Final (9 1/2 min): /*= -6° C (21 °F); n = .58 

Full (10 min): /* = -6°C (21°F); n = . 58 

(a) Tests with V— 67 m/s (150 mph ), MVD = 20 pm 
and I.JFC = .5 g/m 3 . 


Initial (1/2 min): V— 67 m/s (150 mph); n— .58 

Final (9 1/2 min): V- 89 m/s (200 mph); n = .40 

Full (10 min): V— 89 m/s (200 mph); n = .40 


(b) Tests with = -6°C (21°F), A/FD = 20 pm 
and LIFC = .5 g/m 3 . 


Figure 11. Effect of Initial Spray on Final Ice Shape for Airfoil. NACA 0012 Airfoil With Chord, 53.3 cm (21 in). 
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